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IMAGE PROCESSING METHOD, APPARATUS, AND STORAGE MEDIUM 
FOR RECOGNITION OF IRRADIATION AREA 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

The present invention relates to an image 
processing method, an image processing apparatus, and a 
storage medium using secondary differences for 
recognition of an irradiation field. 

10 Related Background Art 

The recent progress in the digital technology 
makes it common to convert a radiographic image to a 
digital image signal, subject this digital image signal 
to image processing, and display it on a CRT or the 

15 like or print it out. In photography of the 

radiographic image it is common practice to diaphragm 
the irradiation area so as to radiate the radiation 
only into a necessary area of a human body on the 
humane grounds or for the reason of preventing 

20 scattering from unnecessary areas so as to prevent 

lowering of contrast. Normally, a processing parameter 
is determined from distribution of density values of 
the image prior to the image processing and the image 
processing is then carried out based on the parameter 

25 thus determined. However, if the irradiation area is 
not limited, unwanted information, in a sense, outside 
an area of interest is used for the determination of 



the image processing parameter, which poses a problem 
that appropriate image processing is not effected. 

It is thus necessary to extract the irradiation 
area and determine the image processing parameter from 
5 information of only the area of interest. As an 

extraction method of the irradiation area, for example, 
Japanese Patent Application (Laid-Open) No. 05-007579 
describes a method for dividing the image area into 
small areas and extracting the irradiation area, based 

10 on values of variance in the small areas. Further, for 
example, Japanese Patent Publication No. 6-90412 
describes a method for expressing a change of image 
densities between a predetermined number of pixels near 
an end portion of the image area by an approximate 

15 equation substantially of a linear equation and 
extracting the irradiation area, based on the 
difference between assumed image density values 
according to this approximate equation and actual 
density values indicated by a sample image signal . 

20 These methods are predicated on the premise that 

an object area in the photographic image includes the 
irradiation area (an area actually irradiated), and a 
process of judging whether the object area is an area 
including the irradiation area (which will also be 

25 referred to as an "area with an irradiation diaphragm" ) 
or an area not including the irradiation area (which 
will also referred to as an "area without the 



irradiation diaphragm" ) is carried out as a preliminary 
operation before execution of these methods. 

Methods for judging presence/absence of the 
irradiation diaphragm include methods for comparing an 
5 average, a median, or the like of densities in a 

central area of the image with an average of densities 
of the object area and judging that the object area is 
an area with the irradiation diaphragm, if the average 
of densities of the object area is not more than a 
10 predetermined value, as described in U.S. Patent No. 
5,091,970, for example. 

SUMMARY OF THE INVENTION 

The method of above Japanese Patent Application 

15 (Laid-open) No. 05-007579 had to obtain the density 

variance in each small area and thus had the problem of 
increased computation complexity and computation time. 
Even within the area of interest, density values vary 
suddenly, for example, at the edge of the lang field 

20 and a change rate of densities could become higher 
there than at the edge of the irradiation area. 
Particularly, the variances become high in the area 
where the lang is in contact with the ribs at the edge 
of the lang. This caused candidate points for the 

25 irradiation area to be extracted from outside the edge 
of the irradiation area as well and this posed another 
problem that the judgment became difficult. 



In the method of , above Japanese Patent Publication 
No. 6-90412, a foot area outside the irradiation area 
was assumed for calculation of the linear approximate 
equation and there arose the problem that the method 
5 was invalid if this foot area was not able to be 

extracted well. Further, it was also assumed that the 
foot area outside the irradiation area included two 
areas of a gentle portion and a quickly rising portion 
and this posed another problem that the method was 
10 invalid where the entire foot area was gentle or where 
there were three or more changes of density gradient. 

The present invention has been accomplished in 
order to solve the problems described above and an 
object of the invention is, therefore, to extract the 
15 irradiation area accurately. 

An image processing method of the present 
invention is an image processing method comprising: 

a step of determining a plurality of areas 
arranged in a predetermined direction on an image and 
20 each having a predetermined shape; 

a step of calculating a secondary difference value 
of density values representing the respective areas in 
the plurality of areas; and 

a step of judging one end point of an irradiation 
25 area from the secondary difference values calculated in 
the calculating step. 

In the conventional image judgment method as 
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described in above U.S. Patent No. 5,091,970 etc., for 
example, where the photographic image is one in a state 
in which the subject overlaps with the edge of the 
image, the average density value at the end of the 
5 image varies depending upon the area of the subject 

over the edge of the image and upon the transmittance 
of the radiation. Therefore, the presence/absence of 
the irradiation diaphragm was sometimes misjudged in 
the area at the end of the image. 

10 When the radiant dose of the irradiation diaphragm 

was small, the density difference also became small 
between the central part of the image and the object 
area and the presence/ absence of the irradiation 
diaphragm in the object area was misjudged in some 

15 cases. 

The present invention has been accomplished to 
eliminate the drawbacks described above and an object 
of the invention is, therefore, to make it possible to 
accurately judge whether an image area as a processed 

20 object is an area including the irradiation area. 

An image processing method of the present 
invention is an image processing method for extracting 
an irradiation area in an input image, the image 
processing method comprising: 

25 a step of detecting an irradiation end, based on a 

density distribution in each area, for a plurality of 
areas in a desired direction in the image; and 
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a step of evaluating the result of the detection, 
based on the result of irradiation ends detected for 
each of the plurality of areas. 

Another image processing method of the present 
5 invention is an image processing method for judging 
whether an object area in an image includes an 
irradiation area, the method comprising; 

a secondary difference value acquisition step of 
acquiring secondary difference values from one- 
10 dimensional image data of the object area; 

an irradiation end extraction step of extracting a 
coordinate of an end point of the irradiation area from 
the secondary difference values acquired in the 
secondary difference value acquisition step; 
15 a comparison step of comparing the coordinate 

extracted in the irradiation end extraction step with a 
coordinate of an end point of the irradiation area 
included in the image, the coordinate being obtained 
preliminarily; and 
20 a judgment step of judging whether the object area 

includes the irradiation area, based on the result of 
the comparison in the comparison step. 

The other objects, arrangements, and effects of 
the present invention will become more apparent with 
25 the description of the embodiments of the present 
invention which follows. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of an irradiation area 
extraction device according to Embodiment 1; 

Fig. 2 is a block diagram to show the structure 
5 inside an irradiation area extraction unit; 

Fig. 3 is a flow chart to show processing 
procedures in the irradiation area extraction unit; 

Fig. 4 is a structural diagram to show calculation 
areas ; 

10 Fig. 5A and Fig. 5B are structural diagrams to 

show a radiographic image and density values of one 
line on the radiographic image; 

Fig. 6A and Fig. 6B are characteristic diagrams to 
show density values of one line on Fig. 5A and values 
15 of secondary difference on this line; 

Fig. 7 is a block diagram to show the structure of 
an image judgment device of Embodiment 2-1; 

Fig. 8 is a diagram for explaining an example of 
image data inputted into the above image judgment 
20 device; 

Fig. 9 is a flow chart for explaining a processing 
program carried out in the above image judgment device; 

Fig. 10 is a block diagram to show the structure 
of an image judgment device of Embodiment 2-2; 
25 Fig. 11 is a flow chart for explaining a 

processing program carried out in the above image 
j udgment devi ce ; 



Fig. 12 is a block diagram to show the structure 
of an image judgment device of Embodiment 3-1; 

Fig. 13 is a diagram for explaining an example of 
image data inputted into the above image judgment 
5 device ; 

Fig. 14 is a flow chart for explaining a 
processing program carried out in the above image 
Judgment device; 

Fig. 15 is a block diagram to show the structure 
10 of an angle extraction device and an area extraction 
device in Embodiment 4-1; 

Fig. 16 is a flow chart of "a processing procedure 
sequence in the angle extraction device and the area 
extraction device; 
15 Fig. 17 is a diagram for explaining the processing 

in the angle extraction device and the area extraction 
device; 

Fig, 18 is a flow chart of a processing procedure 
sequence for obtaining an approximate center in an area 
20 in Embodiment 4-4; 

Fig. 19 is a block diagram to show the structure 
of an area extraction device in Embodiment 4-7; 

Fig. 20 is a flow chart of a processing procedure 
sequence in the area extraction device of Embodiment 4- 
25 7; 

Fig. 21 is a block diagram to show the structure 
of an area extraction device in Embodiment 4-9; and 
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Fig. 22 is a flow chart of a processing procedure 
sequence in the area extraction device of Embodiment 4- 
9. 



5 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

( Embodiment 1 ) 

Fig. 1 is a block diagram to show the structure of 
an irradiation area extraction device according to an 
embodiment of the present invention. In the figure, 

10 reference numeral 101 designates an image input unit 
for accepting input of an image from the outside, 102 
an irradiation area extraction unit for extracting an 
irradiation area irradiated with radiation from the 
image inputted at the image input unit 101, 103 an 

15 image processing unit for performing image processing 
based on the irradiation area extracted at the 
irradiation area extraction unit 102, and 104 a display 
unit for displaying an image subjected to the 
processing at the image processing unit 103 . 

20 Fig. 2 is a block diagram to show the inside 

structure of the irradiation area extraction unit 102, 
in which reference numeral 200 designates a calculation 
area input unit for accepting input of a direction, a 
start point, and an end point for determination of 

25 calculation areas and 201 a calculation area 

determination unit for determining, based on the input 
information at the calculation area input unit 200, 
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areas of calculation carried out by a calculation unit 
202 from the input image supplied from the image input 
unit 101. Numeral 202 denotes a calculation unit for 
calculating primary difference values and a secondary 
5 difference value, described hereinafter, from the 

calculation areas determined at the calculation area 
determination unit 201, 204 a judgment unit for judging 
irradiation area ends, 203 a memory unit for storing 
the values calculated at the calculation unit 202 and 

10 the irradiation area ends judged at the judgment unit 

204, and 205 an irradiation area determination unit for 
determining the irradiation area from the irradiation 
area ends stored in the memory unit 203 and Judged at 
the judgment unit 204. 

15 Fig. 3 is a flow chart to show the flow of the 

processing at the irradiation area extraction unit 102. 
Fig. 4 is a diagram to show the calculation areas 
determined at the calculation area determination unit 
201. Fig. 5A is a diagram to show a radiographic image 

20 and Fig. 5B is a diagram to show density values on line 
Xq to X3 of Fig. 5A. The abscissa indicates positions 
on the line Xq to X3 and the ordinate density values on 
the line. Fig. 6A is an enlarged view of part of Fig. 
5B. Fig. 6B is a plot of secondary difference values 

25 calculated at the calculation unit 202 against points 
on the above line between Xg and X3. 

The operation will be described below. 



- 11 - 



In Fig. 1, an image from an external device not 
illustrated is inputted through the image input unit 
101 into the irradiation area extraction unit 102. The 
irradiation area extraction unit 102 extracts the 
5 irradiation area from the input image and delivers the 
information of the irradiation area to the image 
processing unit 103. The image processing unit 103 
performs the image processing based on the information 
of the irradiation area extracted at the irradiation 

10 area extraction unit 102 and the image through the 

image processing is displayed on the display unit 104. 

The processing at the irradiation area extraction 
unit 102 of Fig. 2 will be described next according to 
the flow of Fig. 3. The irradiation area extraction 

15 unit 102 receives the input of the calculation 
direction, the calculation start point, and the 
calculation end point necessary for the determination 
of the calculation areas at the calculation area input 
unit 200 (step S301 ) . It is, however, noted that the 

20 input herein is required only for an initial input 
screen and it does not have to be set for the next 
input screen and after. Next, the calculation area 
determination unit 201 determines the calculation 
areas, based on the information from the calculation 

25 area input unit 200. The calculation areas herein are 
three areas A, B, and C of a rectangular shape arranged 
in parallel on the image, as illustrated in Fig. 4, and 
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are used as calculation areas for calcula-tion of the 
secondary difference value at a calculation point 
indicated by (x,y). These three areas A, B, and C have 
an equal area which is determined by two parameters of 
5 "a" and "b" and the distance between the areas is 

represented by "d" . These parameters "a", "b", "d" are 
determined experimentally. The calculation area 
determination unit 201 preliminarily determines the 
calculation areas for all calculation points and the 

10 result, together with information of duplicate 

calculation areas, is stored in the memory unit 203 
(step S302). 

Then the calculation unit 202 calculates 
representatives of density values in the above 

15 calculation areas A, B, and C according to a method 
described hereinafter and further calculates the 
primary difference values and the secondary difference 
value. The primary difference values are two values, 
which are a value "e" resulting from subtraction of the 

20 density representative in the area B from the density 
representative in the area C and a value "f" resulting 
from subtraction of the density representative in the 
area A from the density representative in the area B. 
The secondary difference value is a value resulting 

25 from subtraction of the value "f" from the value "e". 

The primary difference values and secondary difference 
value calculated in this way are stored in the memory 
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uni-t 203 . From the duplicate information of the 
calculation areas stored in the memory unit 203, the 
primary difference values for duplicate calculation 
areas are obtained from the values stored in the memory 
5 unit 203 and the secondary difference value thereof is 
calculated using them (steps S303 and S304). 

Then the judgment unit 204 judges the irradiation 
area ends from the primary difference values and the 
secondary difference values stored in the memory unit 

10 203. Fig. 6B shows the secondary difference values 

calculated at the calculation unit 202. The secondary 
difference values are negatively large at the points 
and X2 where the densities vary suddenly, and the 
aforementioned primary difference values f take 

15 positive values in a density increasing direction (when 
seen from the side of the point X^) but take negative 
values in a density decreasing direction. From this 
property the judgment unit 204 judges that a point 
where the above secondary difference value is minimum 

20 and the primary difference values are positive is a 

candidate for an end point of the irradiation area. If 
there are plural candidates, the judgment unit 
determines that a first appearing candidate is an end 
point of the irradiation area ( step S305 ) . Then the 

25 irradiation area end point is stored in the memory unit 
203 (step S306). 

Further, the judgment unit 204 judges whether the 
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processing is completed for all the directions supplied 
from the calculation area input unit 200 and, if not, 
the processing is repeated from step S302 (step S307). 
After the irradiation area end points are determined in 
5 all the directions, the irradiation area determination 
unit 205 determines the irradiation area. Here, a line 
passing the end point of the irradiation area and being 
perpendicular to the calculation direction of the 
secondary difference value (for example, the direction 

10 from Xg to X3 in Figs. 5A and 5B) is calculated for all 
the end points of the irradiation area and an area 
surrounded by these lines obtained is determined as an 
irradiation area (step S308). 

The present embodiment presents the effects of 

15 less computational complexity and capability of readily 
extracting the end point of the irradiation area, 
because the calculation unit 202 is arranged to 
calculate the representatives of density values in the 
calculation areas A, B, C and calculate the secondary 

20 difference value based on the representative values. 

Since the device of the present embodiment has the 
irradiation area determination unit 205 for determining 
the irradiation area, the present embodiment has such 
an effect that the irradiation area can be extracted 

25 even if the irradiation area is rectangular, polygonal, 
or circular. The present embodiment further has such 
an effect that when the irradiation area is polygonal. 
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an almost circumscribed area can be ex-tract:ed by 
selecting directions of lower dimensions than it. 

Described below are methods for calculating the 
representative value in each of the above calculation 
5 areas A, B, C. Simple methods for such calculation 
include a method for calculating an average of all 
density values (all pixel values) in a calculation 
area, a method for sorting all the density values in a 
calculation area according to the density values and 

10 using a density value of a pixel at an intermediate 
position, i.e., a median, and so on. There are also 
methods for using the limited number of pixels in each 
calculation area, instead of all the pixels, by 
sampling or the like according to the necessity and 

15 using an average of density values of the limited 

number of pixels or a median of the limited number of 
pixels . 

Here, when the method for calculating the 
representative of density values in the calculation 

20 area is the method using the average of all the density 
values in the calculation area, the method has the 
effect of being resistant to noise. When the 
representative is determined by the method for sorting 
all the density values in the calculation area 

25 according to the density values and using the density 
value of the pixel at the intermediate position, the 
method has the effects of simple calculation and being 
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resistant to noise. Further, when the representative 
is determined by the method using the average of the 
limited number of points in the calculation area or the 
method for sorting the density values of the limited 
5 number of points and using the density value at the 
intermediate position as a representative value, the 
methods have the effect of capability of further 
decreasing the computational complexity. 

Another approach for determining a representative 

10 of each calculation area A, B, C shown in Fig. 4 will 
be described below. First, let us assume that there 
are the calculation areas having the vertical height of 
"a" and the horizontal height of "b", the distance 
between the areas is represented by "d", and each 

15 calculation point is indicated by (x,y), as illustrated 
in Fig. 4. 

In the present embodiment, where the 
representative value of each area A, B, or C is S(A), 
S(B), or S(C), respectively, the secondary difference 

20 value SS(X) is calculated according to Eq. (1) below. 
SS(X) = S(A) - 2 X S(B) + S(C) ... (1) 

Next described is a method for calculating the 
representative value according to this approach where 
the representative of each area is indicated by S(X) in 

25 general. First, a pixel value at each calculation 

point (x,y) is defined as f(x,y) and values Fl(x) are 
calculated according to Eq. (2) below. 
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Fl (x) = f*f (x,y)dy ' * ' 

Jo 

Here, the integration range is the vertical height 
"a" of the above rectangular region and for simplicity 
of the description thereof Eq. (2) is expressed so as 
5 to perform the integration in the range of "0" to "a". 
Further, values F2(x) are next calculated 
according to Eq. (3) below. 

F2(x) = min{Fl(x + xl ) - h(xl)|xl6K} 
... (3) 

10 Then S(x) is calculated using the values F2(x), 

according to Eq. (4) below. 

S(X) = max{F2(x - xl ) + h(xl)|xieK} 
... (4) 

Here, h(x) is a function represented by Eq. (5) 
15 below and K is a domain of definition thereof. 
h(x) = 0 , -b/4 < X < b/4 

= -00 , otherwise ... ( 5 ) 

Then each of the left and right irradiation ends 
is obtained from a value of x taking minSS(X) in the 
20 left half and the right half of the image. The upper 
and lower ends can also be obtained in a similar 
fashion- 
In other words, the above approach is nothing but 
an operation in which projection of pixel values is 
25 made with respect to a > predetermined direction (the 

vertical direction) of the area, values obtained by the 
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projection are smoothed using a one-dimensional 
morphology filter, and the secondary difference value 
is calculated with the distance d. 

With this calculation method, the irradiation area 
5 can be detected with less errors than in the case of 
the above methods using the average and median. 

The functional blocks 101 to 104, 200 to 205 in 
Fig. 1 and Fig. 2 may be configured on a hardware basis 
or may be configured in a microcomputer system 

10 comprised of CPU, memory, and so on. When they are 
configured in the microcomputer system, the above 
memory comprises the memory medium according to the 
present invention and this memory medium stores a 
program for carrying out the processing illustrated in 

15 the flow chart of Fig. 3. This memory medium can be 

selected from semiconductor memories such as ROM, RAM, 
and the like, optical disks, magneto-optical disks, 
magnetic media, and so on, or it may be substantiated 
in the form of either one of a CD-ROM, a floppy disk, a 

20 magnetic tape, a non-volatile memory card, and so on. 
Since the device according to the present 
embodiment is arranged to determine the calculation 
areas comprised of the plural areas of the 
predetermined shape arranged in the predetermined 

25 direction, calculate the secondary difference value of 
the density values representing the respective areas in 
the plural areas, and judge one end point of the 
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irradiation area from the secondary difference value 
thus calculated, the present embodiment presents the 
effects of capability of decreasing the computational 
complexity and in turn decreasing the computation time 
5 and capability of extracting the end points of the 
irradiation area with accuracy. 

Since the device of the present embodiment is 
arranged to determine the irradiation area from the 
plural end points of the irradiation area judged, it 

10 can enjoy the effect of capability of extracting the 
irradiation area even if the irradiation area is 
rectangular, polygonal, or circular. When the 
irradiation area is polygonal, the device presents the 
effect of capability of extracting the almost 

15 circumscribed area by selecting directions of lower 
dimensions than it. 

Since the density value representing each area in 
the plural areas is an average density value in each 
area, the device can enjoy the effect of capability of 

20 extracting the end points of the irradiation area with 
accuracy even from an image with noise. 
( Embodiment 2 - 1 ) 

The present embodiment is carried out, for 
example, by an image judgment device 2100 illustrated 

25 in Fig. 7. This image judgment device 2100 is provided 
with a judgment unit 2110, a control unit 2120 for 
carrying out control of the operation of the judgment 
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unit 2110, and a program memory 2130 to which the 
control unit 2120 makes access, as illustrated in Fig. 
7. 

The judgment unit 2110 is composed of a secondary 
5 difference value calculation circuit 2201 for 

calculating the secondary difference value from data of 
an object area in an input image, a left end point 
extraction circuit 2202 for extracting a left end point 
of an irradiation area included in the object area, 

10 based on the secondary difference value calculated at 
the secondary difference value calculation circuit 
2201, a right end point extraction circuit 2203 for 
extracting a right end point of the irradiation area 
included in the object area, based on the secondary 

15 difference value calculated at the secondary difference 
value calculation circuit 2201, and a diaphragm 
presence/absence judgment circuit 2204 for judging 
whether the object area is an area with an irradiation 
diaphragm or an area without an irradiation diaphragm, 

20 from the left end point extracted at the left end point 
extraction circuit 2202 and the right end point 
extracted at the right end point extraction circuit 
2203. 

Here, for example, an image G as illustrated in 
25 Fig. 8 is supplied as an input image to the judgment 
unit 2110. This input image G is a two-dimensional, 
radiographic, thoracic part image obtained by 
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photographing the thoracic part with the irradiation 
diaphragm by use of an imaging device with the 
irradiation diaphragm function. 

In above Fig. 8, "X^" and "Xj^" represent positions 
5 of end points ( irradiation area ends ) of the 

irradiation area preliminarily obtained with respect to 
the horizontal axis H, and "A" to "D" are areas at the 
edges of the image. 

These areas A to D are object areas for which 

10 whether the irradiation diaphragm is present or absent 
(i.e., whether the irradiation area is present or 
absent) is judged. In this example, the areas A and B 
are areas with the irradiation diaphragm while the 
areas C and D are areas without the irradiation 

15 diaphragm. 

Various processing programs for controlling the 
operation of the judgment unit 2110 are preliminarily 
stored in the program memory 2130. 

Specifically, for example, a processing program 

20 according to the flow chart as illustrated in Fig. 9 is 
preliminarily stored in the program memory 2130 and 
this processing program is read and executed by the 
control unit 2120, so as to operate the judgment unit 
2110 as follows. 

25 First, when the area C is an object area, the 

secondary difference value calculation circuit 2201 
calculates the secondary difference values SS(x) from 
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data of the area C according to Eq. (21) below (step 
S2301) . 

SS(x) = f(x - d) - 2 X f(x) + f(x + d) 
... (21) 

5 In this Eq. (21), "f(x)" indicates image data of a 

one-dimensional line which crosses the area C in the 
horizontal direction, "x" coordinate values thereof, 
and "d" a constant indicating a difference distance. 

Next, using the secondary difference values SS(x) 
10 calculated at the secondary difference value 

calculation circuit 2201, the left end point extraction 
circuit 2202 extracts the left end point xl of the 
irradiation area included in the area C according to 
Eq. (22) below (step S2302). 
15 SS(xl) = min{SS (x) 1 0 5 X s x} ••• (22) 

x: coordinates on the horizontal axis between the 
left and right ends of the object area (coordinates 
between xa and xb) 

On the other hand, using the secondary difference 
20 values SS(x) calculated at the secondary difference 
value calculation circuit 2201, the right end point 
extraction circuit 2203 extracts the right end point x2 
of the irradiation area included in the area C 
according to Eq. (23) below (step S2303 ) . 
25 SS(x2) = inin{SS(x) tx i x s Length} 

... (23) 

In this Eq. (23), "Length" indicates the length of 
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•the area C along the horizontal axis- 

Then the diaphragm presence/absence judgment 
circuit 2204 obtains orders of approximation of the 
left end point xl extracted at the left end point 
5 extraction circuit 2202 and the right end point x2 
extracted at the right end point extraction circuit 
2203 with respect to the left and right irradiation end 
points preliminarily obtained (the positions and X;, 
of the irradiation ends along the horizontal axis H) 

10 and, if the orders of approximation are high, the 

diaphragm presence/ absence judgment circuit 2204 judges 
that the left end point xl and the right end point x2 
are the left and right end points of the irradiation 
area (the area irradiated directly), i.e., that the 

15 area C is an area without the irradiation diaphragm. 

On the other hand, if the orders of approximation are 
low, it is judged that the area C is an area with the 
irradiation diaphragm (step S2304). 

Specifically, using a constant "e" indicating an 

20 approximation width, if the below conditions are 

satisfied the area C is judged as an area without the 
irradiation diaphragm; otherwise the area C is judged 
as an area with the irradiation diaphragm, 
-e < xl - < e and 

25 -e < x2 - X^ < e 

When in an object area the irradiation ends are 
detected similarly to the irradiation ends in the 
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adjacent area, the object area is judged as an area 
irradiated with radiation as the adjacent area was. 

Steps S2301 to S2304 are carried out for each of 
the other areas A, B, D, in a similar fashion to the 
5 area C. 

As described above, the present embodiment is 
arranged to calculate the secondary difference values 
SS(x) from the data of the object area subjected to the 
judgment of presence/absence of the radiation diaphragm 

10 in the two-dimensional input image G and judge whether 
the object area is an area with or without the 
irradiation diaphragm, using the secondary difference 
values SS(x). 

Construction of the device in this structure 

15 permits boundary points between the directly irradiated 
area and the other areas to be extracted with accuracy 
even from photographic image obtained when a subject 
with low transmittances of radiation is photographed. 
Therefore, presence or absence of the irradiation 

20 diaphragm can be judged with accuracy for the object 

areas including the irradiation area. The presence or 
absence of the irradiation diaphragm can also be judged 
for an object area in a photographic image in which a 
portion with low radiation transmittances, such as the 

25 abdominal part or the like, overlaps with the edge of 
the image - 

It may also be contemplated that on the occasion 
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of extracting the left end point xl and the right end 
point x2 in the object area at the left end point 
extraction circuit 2202 and at the right end point 
extraction circuit 2203, another condition for the 
5 secondary difference values SS(x) used in the 

extraction, for example such a condition that the 
secondary difference values SS(x) are not more than a 
fixed threshold, is added. 
(Embodiment 2-2) 
10 The present embodiment is arranged to apply, for 

example, the projection in the object area represented 
by Eq. (24) below to Eq. (21) in above Embodiment 2-1. 

f{x) = JJf(x,y)dy ^^^^ 

In this Eq. (24), "f(x,y)" represents the image 
15 data of the object area in the input image G, and "x" 
and "y" are coordinates on the horizontal axis and on 
the vertical axis, respectively. Further, "b" and "c" 
indicate the object area. 

When the projection in the object area represented 
20 by above Eq. (24) is used in above Eq. (21) in 

Embodiment 2-1, it can achieve the effects similar to 
those in the case where steps S2301 to S2304 described 
above are carried out in the averaged state of the data 
of the object area. Namely, mis judgment of presence or 
25 absence of the irradiation diaphragm can be prevented 
from being caused by alteration of coordinates of the 
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objec-t area due "to influence of scattered rays, noise, 
and so on. Therefore, presence/ absence of the 
irradiation diaphragm in the object area can be judged 
with better accuracy. 
5 (Embodiment 2-3) 

The present embodiment is applied, for example, to 
an image judgment device 2400 as illustrated in Fig. 
10. 

This image judgment device 2400 has the structure 

10 similar to that of the image judgment device 2100 in 

Embodiment 2-1 described above, but is different in the 
internal structure of the judgment unit 2110 therefrom. 

In the image judgment device 2400 of Fig. 10, 
elements operating similarly to those of the image 

15 judgment device 2100 of Fig. 7 are denoted by the same 
reference symbols and the detailed description thereof 
is omitted herein. Only the different structure from 
Embodiment 2-1 will be described in detail herein. 
The judgment unit 2140 is composed of, as 

20 illustrated in Fig. 10, a coordinate indication circuit 
2404, a secondary difference value calculation circuit 
2401 for calculating the secondary difference values 
from data of an object area (either one of the areas A 
to D in the input image G of Fig. 8 in this example) in 

25 the input image according to coordinates indicated by 
the coordinate indication circuit 2404, a left end 
point extraction circuit 2402 for extracting a left end 
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point of the irradiation area included in the object 
area, based on the secondary difference values 
calculated at the secondary difference value 
calculation circuit 2401, and a right end point 
5 extraction circuit 2403 for extracting a right end 

point of the irradiation area included in the object 
area, based on the secondary difference values 
calculated at the secondary difference value 
calculation circuit 2401. 

10 The judgment unit 2140 further has a memory 

circuit 2405 for storing coordinates of a plurality of 
left and right end points extracted for a plurality of 
lines by the left end point extraction circuit 2402 and 
the right end point extraction circuit 2403, an average 

15 coordinate calculation circuit 2405 for calculating an 
average coordinate at each end of the coordinates of 
the left and right end points stored in the memory 
circuit 2405, a first judgment circuit 2407 for judging 
whether the object area is an area with the irradiation 

20 diaphragm or an area without the irradiation area, 

using the averages of the coordinates of the left and 
right end points calculated at the average coordinate 
calculation circuit 2406, a coordinate variance 
calculation circuit 2408 for calculating variances of 

25 the coordinates of the left and right end points 

according to the judgment result of the first judgment 
circuit 2407, and a second judgment circuit 2409 for 
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again judging whether the object area is an area with 
the irradiation diaphragm or an area without the 
irradiation diaphragm, based on the variances 
calculated at the coordinate variance calculation 
5 circuit 2408. 

The coordinate indication circuit 2404 is 
configured to give an indication of a coordinate for 
calculation of the secondary difference value at the 
secondary difference value calculation circuit 2401 

10 thereto after the left end point extraction circuit 
2402 and the right end point extraction circuit 2403 
extract the left and right end points. 

A processing program, for example, according to 
the flow chart as illustrated in Fig. 11, for 

15 controlling the operation of the judgment unit 2140 is 
preliminarily stored in the program memory 2130 and 
this processing program is read and executed by the 
control unit 2120, whereupon the judgment unit 2140 
operates as follows. 

20 First, supposing the area C is an object area, the 

secondary difference value calculation circuit 2401 
calculates the secondary difference values SSi(x) from 
the data of the area C according to Eq. (25) below for 
each coordinate i indicated by the coordinate 

25 indication circuit 2404 (step S2501 ) . 

SSi(x) = fi(x - d) - 2 X fi(x) + fi(x + d) 
... (25) 
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In this Eq. (25), "fi(x)" represents image data of 
a one-dimensional line crossing the area C in the 
horizontal direction, "x" coordinates thereof, and "i" 
a coordinate of the horizontal line (a coordinate of 
5 the one-dimensional data) indicated by the coordinate 
indication circuit 2404. Further, "d" is a constant 
indicating a difference distance. 

Using the secondary difference values SSi(x) 
calculated at the secondary difference value 
10 calculation circuit 2401, the left end point extraction 
circuit 2402 then extracts the left end point xLi 
according to Eg. (26) below (step S2502). 
SSi (xLi) = minCSSi (x) t 0 < x < x} 
... (26) 

15 x: coordinates on the horizontal axis between the 

left and right ends of the object area (coordinates 
between xa and xb) 

Then the memory circuit 2405 stores the left end 
point xLi extracted at the left end point extraction 
20 circuit 2402 (step S2503). 

Using the secondary difference values SSi(x) 
calculated at the secondary difference value 
calculation circuit 2401, the right end point 
extraction circuit 2403 also extracts the right end 
25 point xRi according to Eg. (27) below (step S2504). 
SSi (xRi) = min{SSi(x) Ix ^ x ^ Length} 
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In this Eq. (27), "Length" indicates the length of 
the area C along the horizontal axis. 

Then the memory circuit 2405 stores the right end 
point xRi extracted at the right end point extraction 
5 circuit 2403 (step S2505). 

After completion of the storage of the left end 
point xLi and the right end point xRi in the memory 
circuit 2405, the coordinate indication circuit 2404 
gives another instruction of a coordinate i of a new 
10 horizontal line to the secondary difference value 
calculation circuit 2401 (step S2506). 

According to this indication, the processing from 
step S2501 is carried out again. This loop processing 
is executed before the coordinate indication circuit 
15 2404 gives an indication of the end to the secondary 
difference value calculation circuit 2401. 

After steps S2501 to S2506 are carried out 
repeatedly and the processing is terminated according 
to the end indication from the coordinate indication 
20 circuit 2404, the memory circuit 2405 is in a storage 
state of the left end points xLi and the right end 
points xRi corresponding to the coordinates i indicated 
during the processing by the coordinate indication 
circuit 2404 to the secondary difference value 
25 calculation circuit 2401. 

Then the average coordinate calculation circuit 
2406 calculates the averages of the left end points xLi 
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and the right end points xRi stored in the memory 
circuit 2405 (step S2507 ) . 

This obtains the average of the left end points 
xLi (left end point average) XL and the average of the 
5 right end points xRi (right end point average) XR. 

Next, similar to the irradiation diaphragm 
judgment circuit 2204 of Fig. 7, the first judgment 
circuit 2407 compares the left end point average XL and 
the right end point average XR calculated at the 

10 average coordinate calculation circuit 2406 with the 
left and right irradiation end points preliminarily 
obtained ( the positions X^ and X^^ of the irradiation 
ends on the horizontal axis H), and if the following 
relations are satisfied it judges that the left and 

15 right end points indicated by the left end point 

average XL and the right end point average XR are the 
left and right end points of the irradiation area, 
i.e., that the area C is an area with the irradiation 
area (or an area without the irradiation diaphragm ). 

20 -e < XL - X^ < e and 

-e < XR - Xj, < e 
On the other hand, if either one is not met, the 
judgment circuit judges that the area C is an area 
without the irradiation area (or an area with the 

25 irradiation diaphragm) (step S2508 ) . 

The next steps S2509 and S2510 are carried out 
only if the first judgment circuit 2407 judges that the 
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area is an area without the irradiation diaphragm. 

Namely, the coordinate variance calculation 
circuit 2408 calculates the variances of the 
coordinates of the left and right end points, using the 
5 left end points xLi and the right end points xRi stored 
in the memory circuit 2405 and the left end point 
average XL and the right end point average XR 
calculated at the average coordinate calculation 
circuit 2406 (step 52509). 
10 This operation provides the variance VL of the 

coordinates of the left end points and the variance VR 
of the coordinates of the right end points. 

The second judgment circuit 2409 compares the 
variances VL and VR of the coordinates of the left and 
15 right end points obtained at the coordinate variance 

calculation circuit 2408 with a predetermined threshold 
TH. When the following relations are satisfied, the 
second judgment circuit 2409 judges that the area C is 
an area with the irradiation area (or an area without 
20 the irradiation diaphragm). 

VL < TH and 

VR < TH 

On the other hand, if either one is not met, the second 
judgment circuit 2409 judges that the area C is an area 
25 without the irradiation area (or an area with the 
irradiation diaphragm) (step S2510)- 

Steps S2501 to S2510 are also carried out for each 
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of the other areas A, B, D in the similar fashion to 
the area C. 

As described above, the present embodiment is 
arranged to carry out the detection of the irradiation 
5 ends on the plural lines in the object area while the 
coordinate indication circuit 2404 gives the 
indications of the coordinates i of one-dimensional 
data lines crossing the object area in the horizontal 
direction to the secondary difference value calculation 

10 circuit 2401. 

Since this structure is arranged to carry out the 
detection of the irradiation ends on the plural lines 
in the object area, the present embodiment is more 
unlikely to be affected by the scattered rays and the 

15 noise and can perform the judgment with better accuracy 
than in the case of the detection of the irradiation 
ends on a certain line. 

Since the variances of the irradiation ends are 
also used for the judgment of presence/ absence of the 

20 irradiation diaphragm by provision of the second 

judgment circuit 2409, for example, when the image end 
portions are exposed to radiation because of the 
influence of the scattered rays, so as to result in 
making the irradiation ends extracted in the image end 

25 portions coincident with the irradiation ends 

preliminarily obtained, the degree of spread of the 
scattered rays can be quantified and thus mis judgment 
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of presence/absence of the irradiation diaphragm can be 
prevented for sure. 
( Embod iment 2 - 4 ) 

The present embodiment is arranged to calculate 
5 the secondary difference values SS(x) according to 

above Eq. (21) and to extract the left end point xl of 
the irradiation area included in the object area, using 
the secondary difference values SS(x), for example, 
similar to Embodiment 2-1 described above. At this 
10 time the present embodiment also uses the sign of the 

primary difference value S(xl) given by Eq. (28) below. 
S(xl) = f(xl) - f(xl - d) ... (28) 

Specifically, for example, if the sign of the 
primary difference value S(xl) is "negative" and if 
15 above Eq. (22) is satisfied, the left end point xl is 
regarded as a left end point of the irradiation area. 

For extracting the right end point x2 of the 
irradiation area included in the object area, the 
present embodiment also uses the sign of the primary 
20 difference value S(x2) given by Eq. (29) below. 

S(x2) = f(x2 + d) - f(x2) ... (29) 

For example, if the sign of the primary difference 
value S(x2) is "negative" and if above Eq. (23) is 
satisfied, the right end point x2 is regarded as a 
25 right end point of the irradiation area. 

Using the left end point xl and the right end 
point x2 obtained in this way, the judgment of 
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presence/ absence of the irradiation diaphragm in the 
object area is carried out in the similar fashion to 
that in Embodiment 2-1 described above. 

When the conditions of the signs of the primary 
5 difference values are also added to the detection of 
the irradiation ends of the object area, 
presence/absence of the irradiation diaphragm in the 
object area can be judged with consideration to the 
inclination of image data outside the irradiation area 

10 due to the scattered rays. This permits the 

presence/absence of the irradiation diaphragm in the 
object area to be judged without mis judgment and with 
better accuracy even if there is a portion or the like 
where densities change quickly in the object area. 

15 (Embodiment 2-5) 

The present embodiment is arranged first to 
calculate the secondary difference values from the data 
of the object area, for example similar to Embodiment 
2-1, but at this time, the present embodiment is 

20 arranged to calculate the secondary difference values 

from data obtained after the data of the object area is 
filtered. 

Specifically, for example, the image data of a 
one-dimensional line of an object area is represented 
25 by "f(x)", the data is subjected to a filtering process 
according to Eqs- (30) and (31) below, and the 
secondary difference values are calculated from the 
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values F2 obtained as a result. 

Fl(x) = min{f(x + xl ) - h(xl)j-d < xl < d} 
... (30) 

F2(x) = max{Fl(x - xl ) + h(xl)|-d < xl < d} 
5 ... (31) 



h(x) { 



= 0 : -d ^ X ^ d 
= -«>: otherwise 



When Embodiment 2-1 is modified so as to calculate 
the second difference values after the data of the 
object area is smoothed by the filtering process as 
described above, the presence/absence of the 
irradiation diaphragm in the object area can be judged 
with better accuracy without influence of the noise, 
particularly, without influence of the noise on the 
line. 

( Embodiment 3 - 1 ) 

The present embodiment is carried out, for 
example, by an image judgment device 3100 as 
illustrated in Fig. 12. 

This image judgment device 3100 is composed of a 
judgment unit 3110, a control unit 3120 for controlling 
the operation of the judgment unit 3110, and a program 
memory 3130 to which the control unit 3120 makes 
access, as illustrated in Fig. 12. 

The judgment unit 3110 is composed of a coordinate 
indication circuit 3201, a characteristic value 
extraction circuit 3202 for extracting a characteristic 
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value from data of an object area in an input image 
according to a coordinate indicated by the coordinate 
indication circuit 3201, an end point extraction 
circuit 3203 for extracting a coordinate of an end 
5 point of the irradiation area included in the object 
area, based on the characteristic value extracted at 
the characteristic value extraction circuit 3202, a 
memory circuit 3204 for storing coordinates of end 
points extracted at the end point extraction circuit 

10 3203, and a judgment circuit 3205 for judging whether 
the object area is an area with the irradiation 
diaphragm or an area without the irradiation diaphragm, 
from the coordinates of the end points stored in the 
memory circuit 3204. 

15 Here, for example, an image G illustrated in Fig. 

13 is supplied as an input image to the judgment unit 
3110. This input image G is a two-dimensional, 
radiographic, thoracic part image obtained by 
photographing the thoracic part with the irradiation 

20 diaphragm by use of an imaging device with the 
irradiation diaphragm function. 

In above Fig. 13, "X^" and "Xj," represent positions 
of irradiation area ends on the horizontal axis X. "A" 
to "D" are areas at the edges of the image part. In 

25 this example, the areas A and B out of these "A" to "D" 
are areas with the irradiation diaphragm while the 
other areas C and D are areas without the irradiation 
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diaphragm. 

Various processing programs for controlling the 
operation of the judgment unit 3110 are preliminarily 
stored in the program memory 3130. 

Specifically, for example, a processing program 
according to the flow chart as illustrated in Fig. 14 
is preliminarily stored in the program memory 3130 and 
this processing program is read and executed by the 
control unit 3120 to operate the judgment unit 3110 as 
follows. 

Described below is a process of judging whether 
the irradiation diaphragm is present or absent at the 
lower edge. 

First, the coordinate indication circuit 3201 
gives indications of plural coordinates on the X-axis 
between X, and X^ in the input image G of above Fig. 13 
to the characteristic value extraction circuit 3202 
(step S3301). 

Specifically, the coordinate indication circuit 
3201 successively indicates, for example, coordinates i 
(= 1 to 10) on the X-axis, which are ten points 
distributed at equal intervals between X^ and X^^. 

According to each coordinate i indicated by the 
coordinate indication circuit 3201, the characteristic 
value extraction circuit 3202 then calculates a 
characteristic value of one-dimensional image data 
corresponding to the coordinate i, for example, the 
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secondary difference values SSi(y) according to Eq. 
(41) below (step S3302 ) . 

SSi(y) - fi(y - d) - 2 x fi(y) + fi(y + d) 
... (41) 

5 In this Eq. (41), "fi(y)" indicates image data of 

a one-dimensional line crossing the object area in the 
vertical direction, "y" coordinates thereof, and "d" a 
constant indicating a difference distance. 

Using the secondary difference values SSi(y) 
10 obtained at the characteristic value extraction circuit 
3202, the end point extraction circuit 3203 then 
extracts a coordinate yi of an end point ( an 
irradiation end point) of the irradiation area included 
in the object area according to Eg. (42) below (step 
15 S3303). 

SS(yi) = min{SSi(y)jO S y S Length} 
... (42) 

In this Eq. (42), "Length" represents the length 
of the input image G along the vertical axis ( the 
20 direction of the Y-axis ) - 

Then the memory circuit 3204 stores the coordinate 
yi of the irradiation end point obtained at the end 
point extraction circuit 3203 (step S3304). 

The above processing of steps S3301 to S3304 is 
25 carried out repeatedly before the coordinate i 

indicated by the coordinate indication circuit 3201 
reaches "10". This results in storing the coordinates 
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yl, y2, y3,..., ylO of the irradiation end points 
corresponding to the coordinates i = 1, 2, 3,..., 10 
indicated by the coordinate indication circuit 3201, in 
the memory circuit 3204. 
5 After completion of the processing of steps S3301 

to S3304 up to the coordinate i = 10, the judgment 
circuit 3205 then calculates an average of the 
coordinates yl, y2, y3,..., ylO stored in the memory 
circuit 3204 (step S3305) and calculates a variance Bv 

10 of the irradiation end points (step S3306 ) . 

Then the judgment circuit 3205 judges that the 
object area is an area without the irradiation 
diaphragm, if the variance Bv is not less than a 
predetermined threshold TH; otherwise, it is judged 

15 that the object area is an area with the irradiation 
diaphragm ( step S3307 ) . 

The threshold TH is a constant which is determined 
experimentally. 

As described above, the present embodiment is 

20 arranged to carry out the extraction of irradiation 
ends of plural column lines in the object area for 
which presence/absence of the irradiation diaphragm is 
judged and judge the presence/absence of the 
irradiation diaphragm of the object area according to 

25 the variance of the irradiation ends. If there is an 
irradiation area in the object area the coordinates of 
the irradiation end points of the plural column lines 
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will be aligned approximately on one horizon-tal axis 
and -the variance thereof will be thus small. If there 
is no irradiation area in the object area the 
coordinates of the irradiation end points of the plural 
5 column lines will be distributed and the variance will 
be thus large. 

Therefore, the presence/absence of the irradiation 
diaphragm in the object area can be determined with 
accuracy because of the structure using the variance. 

10 Since the device is constructed to use the 

secondary difference value in order to extract the 
irradiation end, boundary points can be extracted with 
accuracy between an area irradiated directly and the 
other areas even in a photographic image obtained by 

15 photographing a subject with low transmittances of 
radiation. Therefore, the device of the present 
embodiment can judge the presence/absence of the 
irradiation diaphragm in the object area including the 
irradiation area with accuracy. In addition, the 

20 presence/absence of the irradiation diaphragm in the 
object area can be judged with accuracy even in a 
photographic image in which a portion with low 
radiation transmittances such as the abdominal part or 
the like overlaps with an end portion of the image. 

25 Although the present embodiment is arranged to use 

the secondary difference value in order to make a 
judgment of the irradiation diaphragm, the apparatus of 
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the present invention does not always have to be 
limited to this; for example, where change of density- 
is quick at an irradiation end, the apparatus may also 
be arranged to use the primary difference values or 
5 higher-order difference values. In this case, the 
primary difference values or the higher-order 
difference values are obtained from the object area and 
a first appearing point of a value not less than a 
predetermined threshold is regarded as a candidate for 

10 an irradiation end. 

Although the present embodiment is arranged to use 
the variance in order to make a judgment of the 
irradiation diaphragm, the apparatus of the present 
invention does not always have to be limited to this; 

15 for example, the apparatus may also be arranged to use 
another index that can indicate discrete degrees of 
frequency, such as standard deviation. 
(Embodiment 3-2) 

The present embodiment is arranged to apply the 

20 projection of a constant width in an object area 

according to Eq. (43) below to Eq. (41) in Embodiment 
3-1 described above. 

fi(y) = r fU,y)dx ^^^^ 

Jxi-l 



25 



In this Eq. (43), "fi(y)" represents image data of 
a one-dimensional line in the object area in above Eq. 
(41), and "Xi" and "Xi-1" represent coordinates 
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indicated by the coordinate indication circuit 3201. 

When Embodiment 3-1 is modified so as to apply the 
projection in the object area represented by above Eq. 
(43) to above Eq. (41) as described above, the 
5 modification can enjoy the same effects as in the case 
where aforementioned steps S3301 to S3307 are carried 
out in an averaged state of the data of the object 
area. Namely, this structure can prevent mis judgment 
about the presence/absence of the irradiation diaphragm 

10 from being caused by change of coordinates of the 

object area due to influence of the scattered rays, the 
noise, and so on. The presence/absence of the 
irradiation diaphragm in the object area can be judged 
with better accuracy accordingly. 

15 (Embodiment 3-3) 

The present embodiment is arranged to calculate 
the secondary difference values SSi(y) according to 
above Eq. (41) and extract the irradiation end of the 
object area using the secondary difference values 

20 SSi(y), for example, as Embodiment 3-1 described above 
was. At this time, the present embodiment also uses 
the sign of the primary difference value Si(y) 
expressed by Eq. (44) below. 

Si(y) = fi(y) - fi(y - d) ... (44) 

25 Specifically, for example, if the sign of the 

primary difference value Si(y) is "negative" and if 
above Eq. (42) is satisfied, the point "y" is regarded 
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as an irradiat;ion end. 

Using the irradiation ends obtained in this way, 
whether the irradiation diaphragm is present or absent 
in the object area is judged in the similar fashion to 
5 aforementioned Embodiment 3-1. 

When Embodiment 3-1 is modified to use the sign of 
the primary difference value as well on the occasion of 
detecting the irradiation end in the object area, the 
presence/absence of the irradiation diaphragm in the 

10 object area can be judged, also taking the inclination 
of image data outside the irradiation area due to the 
scattered rays into consideration. Therefore, whether 
the irradiation diaphragm is present or absent in the 
object area can be judged with better accuracy without 

15 misjudgment, even if there is a quickly changing 

portion of density or the like within the object area. 
(Embodiment 3-4) 

The present embodiment is arranged first to 
calculate the secondary difference values from the data 

20 of the object area, similar to Embodiment 3-1; but at 
this time, the data of the object area is subjected to 
a filtering process and the secondary difference values 
are calculated from the data after the filtering 
process. 

25 Specifically, for example, where the image data of 

a one-dimensional line of the object area is "f{x)", it 
is subjected to the filtering process according to Eq. 
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(45) and Eq. (46) below and the secondary difference 
values are calculated from values F2 obtained as a 
result . 

Fli(y) = inin{fi(y + yl ) - h(yl)|-d < yl < d} 
5 ... (45) 

F2i(y) = max{Fli(y - yl ) + h(yl)|-d < yl < d} 
... (46) 

y,(^^ ^ r 0 : "d < X S d 

^^^-^^ i-oo; otherwise 

When Embodiment 3-1 is modified so as to calculate 

10 the secondary difference values after the data of the 
object area is smoothed by the filtering process as 
described above, whether the irradiation diaphragm is 
present or absent in the object area can be judged with 
better accuracy without being affected by the noise, 

15 particularly, without being affected by the noise on 
the line. 
( Embodiment 4-1 ) 

Fig. 15 is a block diagram to show the structure 
of an angle extraction device and an area extraction 

20 device of two-dimensional image data according to 

Embodiment 4-1. In Fig. 15, reference numeral 4101 
designates a characteristic quantity calculation means 
for calculating a characteristic quantity from data of 
a one-dimensional line at a coordinate indicated by a 

25 coordinate indication means 4103, 4102 an end point 
extraction means for extracting an area end portion 
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from -the characteristic quantity calculated at the 
characteristic quantity calculation means 4101, 4103 a 
coordinate indication means for indicating a coordinate 
of image data to be calculated by the characteristic 
5 quantity calculation means 4101, 4104 an end point 

memory means for storing a coordinate of an end point 
extracted at the end point extraction means 4102, and 
4105 a rotation angle indication means for indicating a 
rotation angle of a rotation axis onto which 

10 coordinates of end points stored in the end point 
memory means 4104 are projected. 

Reference numeral 4106 denotes an accumulated 
quantity calculation means for calculating an 
accumulated quantity of the coordinates of the end 

15 points projected onto the aforementioned rotation axis 
and stored in the end point memory means 4104, 4107 an 
accumulated quantity memory means for storing the 
accumulated quantity calculated at the accumulated 
quantity calculation means 4106, 4108 a rotation angle 

20 judgment means for judging a rotation angle from the 

accumulated quantity stored in the accumulated quantity 
memory means 4107, and 4109 an end point judgment means 
for extracting an area end from the angle judged at the 
rotation angle judgment means 4108 and the accumulated 

25 quantity stored in the accumulated quantity memory 
means 4107. 

Fig. 16 is a flow chart of a processing sequence 
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in the angle ex-traction device and the area extraction 
device in Embodiment 4-1. Fig. 17 is a schematic 
diagram for explaining the processing in the angle 
extraction device and area extraction device of the 
5 present embodiment, in which "a" represents two- 
dimensional data (data within a light receiving surface 
in the case of copiers etc., or data in a sensor 
surface in the case of X-ray apparatus etc.), "b" a 
rectangular area (an irradiation area in the case of 

10 the X-ray apparatus or a sheet surface in the case of 
the copiers etc.), "c" a start point of the rotation 
axis, "d" an area end point, "e" a start axis of angle, 
"f" a rotation axis, and "g" an angle difference Gn 
between the start axis e of angle and the rotation axis 

15 f . 

The operation of the present embodiment will be 
described below according to the flow chart of Fig. 16. 

The characteristic quantity calculation means 4101 
of Fig. 15 calculates the secondary difference values 
20 according to Eq, (51) below (steps S4101 and S4102). 

In the equation fi(x) represents image data of the i-th 
row line indicated by the coordinate indication means 
4103, SSi(x) the secondary difference values, and "c" a 
constant . 

25 SSi(x) = fi(x - c) - 2 X fi(x) + fi(x + c) 

... (51) 

Then the end point extraction means 4102 extracts 
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a candidate for an end point: of the area according to 
Eq. (52) below (step S4103). 

For the left area end of the i-th row, the 
following equation is used. 
5 xi = min{SSi (x) I 0 s x < x} ... (52) 

In this equation x indicates a coordinate on an axis 
along an indicated direction within the area (the 
details will be described hereinafter and reference 
should be made to Eq. (76)). Then the end point memory 

10 means 4104 (Fig. 15) stores the coordinate of the area 
end obtained according to Eq. (52) (step S4104), Then 
steps S4101 to S4104 are repeated before the processing 
is complete for all rows indicated by the coordinate 
indication means 4103 (Fig. 15). 

15 Next, the rotation angle indication means 4105 

(Fig. 15) indicates the rotation angle 9n (Fig. 17) 
(step S4105), and the accumulated quantity calculation 
means 4106 (Fig. 15) performs processing indicated by 
Eq. (53) to Eq. (57) below to calculate the accumulated 

20 quantity L(9n,X) of coordinates x on the aforementioned 
rotation axis f (Fig. 17) (step S4106). The projection 
of the area end points onto the aforementioned rotation 
axis is obtained from projected values of candidate 
points onto the rotation axis rotated about the start 

25 point c(xc,yc) of the aforementioned rotation axis. 

Xi,„ = -(Xi - xc)*cos(en) + (Yi - yc)*sin(en) 
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Here, Xi, Yi represent coordinates of a candidate point 
for the area end of the i-th row, and X^^^ represents 
projected positions of the candidate points for the 
area ends on the rotation axis x, which are points 
resulting from rotation by Bn about the start point c 
(Fig. 17) of the aforementioned rotation axis. 

L(en,x) =J^Z(X,Xien) ... (54) 

i=i 

Here, L(en,X) indicates the number of projected 
candidate points on the axis resulting from the 
rotation of the aforementioned rotation axis by en. 
Here, Z{X,Xi) = 1, X - dl < Xi < X + dl 

= 0, otherwise, ... (55) 

where dl is a constant. 

Then the accumulated quantity memory means 4107 
(Fig. 15) stores the accumulated quantity L(en,x) 
obtained by Eq. (54) (step S4107). Further, the 
processing of steps S4105 to S4107 is repeated for all 
angles en indicated by the rotation angle indication 
means 4105 (Fig. 15). 

Then the rotation angle judgment means carries out 
the processing of Eq. (56) to Eq. (58) below to extract 
the rotation angle of the area (step S4108). 

n(en) = L(en,Xmax - 1) + L(en,Xmax) + L(en,Xmax + 1) 
... (56) 

This represents the number of overlaps of area ends at 
the maximum overlap point in the projection of area 
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ends onto the rotation axis. Here, Xmax indicates a 
coordinate of X that satisfies the following. 

L(en,Xmax) = max{L( en, X ) | xekl } ... (57) 

Here, kl represents coordinates on the rotation axis. 
5 After all, the rotation angle 6 is an angle 

obtained according to the following equation. 

G = inax{n( en) ] en£K} ... (58) 

In this equation K is an arbitrary domain of 
definition. 

10 Then the end point Judgment means 4109 (Fig. 15) 

extracts Xmax satisfying Eq. (59) as a candidate for a 
left end point of the rectangular area from 9 
determined by Eq. (58) and the accumulated quantities 
L(e,X) stored in the accumulated quantity memory means 

15 4107 (Step S4109). 

L(e,Xmax) = max{L( e,X) I xekl} ... (59) 
Likewise, candidates are also extracted for the 
right, upper, and lower end points. It is, however, 
noted that they do not always have to be extracted, 

20 because the angle of the area is already determined. 

As described above, the present embodiment has the 
following effect because of the use of the secondary 
difference value; when the area is extracted, boundary 
points between the area directly irradiated with X-rays 

25 and the other areas can be extracted as end points of 
the area with accuracy even from an object with low X- 
ray transmittances . 
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Further, since the angle and end points of the 
area are determined from the accumulation of end 
points, the angle and end points of the area can be 
extracted with accuracy and within short computation 
5 time. 

In the case of the copiers, FAX, OCR, etc., a 
necessary area can be extracted based on the area and 
angle obtained. Further, the present embodiment has 
the effect of capability of automatically carrying out 
10 selection of a sheet with accuracy, because the size of 
the area can be identified. In addition, even if a 
sheet is placed obliquely, correction can be made and 
thus the result is equivalent to that where the sheet 
is placed at a normal position. 
15 (Embodiment 4-2) 

In the present embodiment, where the image data of 
a one-dimensional line is f(x) and the secondary 
difference values thereof are SS(x) defined by Eq. (1), 
the sign of the primary difference value S(XL) of Eq. 
20 (60) below is also added to the operation of extracting 
the end point XL at the end point extraction means 4102 
(Fig. 15). 

S(XL) = f(XL) - f(XL - d) ... (60) 

For example, for extracting the left end point, XL 
25 that makes S(XL) negative and that satisfies Eq. (61) 
below is regarded as a left end point. 

SS(XL) = min{SS (x) I 0 ^ X ^ x} (61) 
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As described above, the present embodiment can 
also take the inclination of the image data outside the 
area into consideration, because the sign of the 
primary difference is added to the extraction of the 
5 area end. This eliminates erroneous extraction of a 
quickly changing portion of density in the area and 
thus the area end point can be extracted with better 
accuracy . 
( Embod iment 4 - 3 ) 
10 In the present embodiment, where the one- 

dimensional data of the i-th row is fi(x) and is 
subjected to the filtering process according to Egs. 
(64) and (65) below and values resulting from the 
filtering process are defined as Fli(x) and F2i(x), the 
15 characteristic quantity calculation means 4101 (Fig. 
15) uses F2i(x) for calculating the secondary 
difference values defined by Eq. (51). 

Fli(x) = min{fi(x + xl ) - h(xl)|-d < xl < d} 
... (64) 

20 Further, F2i(x) is defined as follows. 

F2i(x) = max{Fli(x - xl ) + h(xl)|-d < xl < d} 
... (65) 

Here, h(x) is a function defined below. 

h(x) =0, -d < xl < d 
25 = oo, otherwise ... (66) 

Hence, the secondary difference values are 
calculated as follows. 
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SSi(x) = F2i(x - c) - 2 X F2i(x) + F2i(x + c) 



(67) 



10 



15 



20 



As described above, the present embodiment is 
arranged to smooth the one-dimensional image data for 
calculation of the secondary difference values by the 
filtering process, thereby accomplishing the effect of 
being not affected by the noise, particularly by the 
noise on the line. 
( Embodiment 4 - 4 ) 

The present embodiment is an example in which the 
start point c (Fig. 17) of the aforementioned rotation 
axis is placed at the barycenter of data not less than 
a fixed density. The flow of processing will be 
described according to a flow chart of a procedure 
sequence for obtaining an approximate center in the 
area, shown in Fig. 18. 

First, an accumulation histogram of whole image 
data is produced in step S4201 and upper densities 
above TH % are extracted (steps S4202, S4203, S4204). 
Then the barycenter of the upper densities above TH % 
is regarded as an approximate barycenter of the 
rectangular area. 




... (68) 




(69) 
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Here, f(x,y) represents an image of the upper 
densities above TH % and x stand for the barycenter. 

As described above, since the present embodiment 
is arranged to calculate the barycenter of data above 
5 the fixed density value, the present embodiment can 
extract the approximate barycenter of the area of 
interest with accuracy without being affected by the 
peripheral data, even if the object area is in the 
peripheral part of the two-dimensional image data. 
10 Further, since the start point c (Fig. 17) of the 

aforementioned rotation axis f is placed in the object 
area, degrees of change in accumulated quantities of 
area ends become large against rotation of the rotation 
axis and the present embodiment can thus enjoy the 
15 effect of capability of extracting the rotation angle 
and area with better accuracy. 
( Embodiment 4- 5 ) 

The present embodiment is arranged as follows . 
When the end point extraction means 4102 (Fig. 15) 
20 extracts an end point, a point is not extracted as a 

candidate for an end point if f maxi < Th at the point, 
f maxi being defined below. 

f maxi = max{f i(x) I xEk} ... (70) 
Here, "k" indicates a domain of definition of row data 
25 and Th a fixed threshold value. 

As described above, the present embodiment can 
eliminate end points that can be points outside the 
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object area with high possibilities and thus can 
extract the area and the angle of the area with better 
accuracy . 
( Embodiment 4 - 6 ) 

The present embodiment uses a first order 
differential or a higher order differential as a 
characteristic quantity at the characteristic quantity 
calculation means 4101 (Fig. 15). When Di(x) 
represents differential values of the i-th row, the end 
point extraction means 4102 (Fig. 15) extracts a point 
satisfying the condition of xi > th as a candidate for 
an end point on the occasion of extracting a candidate 
for the left end point of the i-th row. 

As described above, the present embodiment has the 
effect of shorter computation time, because it adopts 
the simple differentiation. 
( Embodiment 4-7 ) 

Fig. 19 is a block diagram to show the structure 
of an area extraction device according to the present 
embodiment. In Fig. 19, reference numeral 4301 
designates a secondary difference value calculation 
means for calculating the secondary difference values 
of one-dimensional image data in a designated direction 
(for example, which is determined according to the 
value determined by the angle extraction device 
described above), 4302 a left end point extraction 
means for extracting a left end point of the area. 
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based on the secondary difference values calculated at 
the secondary difference value calculation means 4301, 
and 4303 a right end point extraction means for 
extracting a right end point of the area, based on the 
secondary difference values calculated at the secondary 
difference value calculation means 4301. Fig. 20 is a 
flow chart of a processing procedure sequence in the 
area extraction device according to the present 
embodiment . 

The flow of the processing in the present 
embodiment will be described according to Fig. 20. The 
secondary difference value calculation means 4301 (Fig. 
19) calculates the secondary difference values SS(x) 
according to a calculation equation defined by Eq. (71) 
below (step S4301). Here, f(x) represents the one- 
dimensional data of a line crossing the area in the 
designated direction and x represents coordinates 
thereof. Further, "d" denotes a constant indicating a 
difference distance. 

SS(d) = f(x - d) - 2 X f(x) + f(x + d) 
... (71) 

The left end point extraction means 4302 (Fig. 19) 
extracts the left end point xl according to Eq. (72) 
below (step S4302 ) . Here, x represents a coordinate on 
the horizontal axis in the area. 

SS (xl) = min{SS (x) I 0 ^ x ^ x} • • • (72) 

Then the right end point extraction means 4203 
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(Fig. 19) extracts the right end point x2 according to 
Eq. (73) below (step S4303 ) . 



SS(x2) = min{SS(x) Ix < x < Length} 



(73) 



Here, "Length" indicates the length of the image data 
along the horizontal axis- 

As described above, since the present embodiment 
uses the secondary difference values, it has the effect 
of capability of extracting the area with accuracy even 
if the density change is gentle at the area end. 
( Embodiment 4 - 8 ) 

The present embodiment uses the projection of the 
image area as defined by Eq. (74) below, as the one- 
dimensional data f(x) in aforementioned Eq. (71). 

f (x) = JJf (x,y)dy ... (74) 




Here, f(x,y) represents the image data in the 
designated direction and x, y coordinates on the 
horizontal and vertical axes, respectively (which are 
inclined in a certain direction). Further, "b", "c" 
represent an arbitrary domain across the object area. 

As described above, since the present embodiment 
uses the projection in the image area, it has the same 
effect as in the case of averaging the data of the 
image area. Therefore, the present embodiment can 
prevent the extraction from being affected by the noise 
inside and outside the object area and thus can extract 
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the area with better accuracy, 
( Embodiment 4- 9 ) 

Fig. 21 is a block diagram to show the structure 
of an area extraction device according to the present 
5 embodiment. In Fig. 21, reference numeral 4401 

designates a secondary difference value calculation 
means for calculating the secondary difference values 
of one-dimensional image data in a designated direction 
( for example, which is determined according to the 

10 value determined by the angle extraction device 

described above) indicated by a coordinate indication 
means 4404, 4402 a left end point extraction means for 
extracting a left end point of the area, based on the 
secondary difference values calculated at the secondary 

15 difference value calculation means 4401, a right end 

point extraction means for extracting a right end point 
of the area, based on the secondary difference values 
calculated at the secondary difference value 
calculation means 4401, and 4404 the coordinate 

20 indication means for indicating a coordinate of one- 
dimensional data for calculation of the secondary 
difference values at the secondary difference value 
calculation means 4401 after extraction of the right 
area end at the right end point extraction means 4403 . 

25 Numeral 4405 denotes a memory means for storing 

coordinates of right end points and left end points 
extracted at the left end point extraction means 4402 
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and at the right end point extraction means 4403, and 
4406 an average coordinate calculation means for 
calculating an average coordinate for each end point 
from the left and right ends stored in the memory means 
5 4405. 

Fig. 22 is a flow chart of a processing procedure 
sequence of the area extraction device according to the 
present embodiment. The operation will be described 
according to Fig. 22. 

10 The secondary difference value calculation means 

4401 (Fig, 21) calculates the secondary difference 
values SSi(x) according to a calculation equation 
defined by Eq. (75) below (step S4401) . Here, fi(x) 
represents one-dimensional data of a line crossing the 

15 area in the designated direction, x coordinates 
thereof, and "i" a coordinate of the row in the 
designated direction indicated by the coordinate 
indication means 4404. Here, "d" represents a constant 
indicating a difference distance. 

20 SSi(x) = fi(x - d) - 2 X fi(x) + fi(x + d) 

... (75) 

Then the left end point extraction means 4402 
(Fig. 21) extracts the left area end XLi according to 
Eq. (76) below (step S4402). 
25 Here, x represents a coordinate on the axis along 

the designated direction in the area. 
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SSi (XLi) = min{SSi (x) I 0 ^ x < x} 
... (76) 

Then the memory means 4405 (Fig. 21) stores the 
left area end XLi extracted at the left end point 
5 extraction means 4402 (step S4403 ) . 

Next, the right end point extraction means (Fig, 
21) extracts the right area end XRi according to Eq. 
(77) below (step S4404 ) - 

SSi (XRi) = min{SSi (x) Ix ^ x ^ Length} 
10 ... (77) 

Here, "Length" represents the length of the image 
data along the horizontal axis in the predetermined 
direction. 

Then the memory means 4405 (Fig. 21) stores the 
15 right area end XRi extracted at the right end point 
extraction means 4403 (step S4404). 

Next, the coordinate indication means 4404 (Fig. 
21) indicates a coordinate of a row in a new designated 
direction and step S4401 to step S4405 are repeated. 
20 An indication of the end of the loop is issued from the 
coordinate indication means 4404 (Fig. 21) (step 
S4406) . 

Next, after completion of the above extraction of 
the left and right end points, the average coordinate 
25 calculation means 4406 (Fig. 21) calculates an average 
of the coordinates of the left and right end points 
stored in the memory means 4405 ( step S4407 ) . 
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As described above, since the present embodiment 
uses the secondary difference values, the boundary 
points of the area can be extracted with accuracy even 
if the density values vary gently at the boundary of 
5 the area. 

Further, since the area end is determined from the 
average of the area ends on plural lines, the present 
embodiment has the effects of being more resistant to 
the noise and capable of extracting the area with 
10 higher accuracy than in the case using a single point. 
(Embodiment 4-10) 

In the present embodiment, where the image data of 
a one-dimensional line in the designated direction is 
expressed by f(x) and the secondary difference values 
15 thereof by SS(x) defined by Eq. (78) below, the sign of 
the primary difference value S(XL) defined by Eg. (79) 
below is added to the extraction of the left end point 
XL at the left end point extraction means 4302 (Fig. 
19 ) or 4402 (Fig. 21) . 
20 SS(x) = f(x - d) - 2 X f{x) + f(x + d) 

... (78) 

S(XL) = f(XL) - f(XL - d) ... (79) 

For example, for extracting a left end point, XL 
that makes s(XL) negative and that satisfies Eq. (80) 
25 below is regarded as a left end point. 

SS (XL) = min{SS (x) f 0 ^ x ^ x} - • • ( 80 ) 
S(XR) = f(XR + d) - f(XR) ... (81) 
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Likewise, where the right end point is extracted 
at the right end point extraction means 4303 (Fig. 19) 
or 4403 (Fig. 21), XR that makes S(XR) defined by above 
Eq. (81) negative and that satisfies Eq. (82) below is 
5 regarded as a right area end. 

SS (XR) = min{SS(x) Ix < x < Length} 
... (82) 

As described above, since the present embodiment 
is arranged to add the sign of the primary difference 

10 to the extraction of the area end, it can also take the 
inclination of the image data outside the area into 
consideration. Therefore, erroneous extraction can be 
prevented at a quickly changing portion of density in 
the area, so that the area can be extracted with better 

15 accuracy. 

( Embodiment 4-11) 

In the present embodiment, where the one- 
dimensional data in the designated direction is defined 
by f(x) and values resulting from the filtering process 

20 according to Eqs . (83), (84) below are defined by 

Fl(x), F2(x), the values F2(x) are used for calculation 
of the secondary difference values at the secondary 
difference value calculation means 4301 (Fig. 19) or 
4401 (Fig. 21). 

25 Fl(x) = min{f(x + xl ) - h(xl)|-d < xl < d} 

... (83) 

Further, F2(x) is defined as follows. 



- 63 - 



F2{x) = max{Fl(x - xl ) + h(xl)l-d < xl < d} 
... (84) 

Here, h(x) is a func-tion defined as follows- 
h{x) =0, -d < xl < d 
5 = -oo^ otherwise ... (85) 

As described above, since the present embodiment 
is arranged to smooth the one-dimensional image data 
for the computation of the secondary difference values 
by the filtering process, it has the effect of being 
10 not affected by the noise, particularly, by the noise 
on the line. 

While the present invention has been described 
with respect to what is presently considered to be the 
preferred embodiments, it is to be understood that the 
15 invention is not limited to the disclosed embodiments. 
To the contrary, the invention is intended to cover 
various modifications and equivalent arrangements 
included within the spirit and scope of the appended 
claims. The scope of the following claims is to be 
20 accorded the broadest interpretation so as to encompass 
all such modifications and equivalent structures and 
functions . 
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WHAT IS CLAIMED IS: 

1. An image processing method comprising: 
a step of determining a plurality of areas 

arranged in a predetermined direction on an image and 
each having a predetermined shape; 

a step of calculating a secondary difference value 
of density values representing the respective areas in 
said plurality of areas; and 

a step of judging one end point of an irradiation 
area from said secondary difference values calculated 
in said calculating step. 

2. A method according to Claim 1, further 
comprising a step of determining said irradiation area 
from a plurality of end points of the irradiation area 
judged in said judging step. 

3. A method according to Claim 1, wherein said 
density values representing the respective areas in 
said plurality of areas are average density values in 
the respective areas, 

4. A method according to Claim 1, wherein said 
density values representing the respective areas in 
said plurality of areas are medians of density values 
in the respective areas. 



- 65 - 



5. A method according to Claim 1, wherein said 
density values representing the respective areas in the 
plurality of areas are averages of density values at a 
limited number of points in the respective areas. 

5 

6. A method according to Claim 1, wherein said 
density values representing the respective areas in the 
plurality of areas are medians of density values at a 
limited number of points in the respective areas. 

10 

7. A method according to Claim 1, wherein said 
density values representing the respective areas in the 
plurality of areas are calculated using integrated 
values in a predetermined direction of pixels in said 

15 plurality of areas. 



8. A method according to Claim 7, wherein said 
density values representing the respective areas in 
said plurality of areas are obtained by smoothing said 

20 integrated values . 

9. An image processing method for extracting an 
irradiation area in an input image, said image 
processing method comprising: 

25 a step of detecting an irradiation end, based on a 

density distribution in each area, for a plurality of 
areas in a desired direction in said image; and 
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a step of evaluating the result of said detection, 
based on the result of irradiation ends detected for 
each of said plurality of areas. 

5 10. A method according to Claim 9, wherein said 

detection is carried out using secondary difference 
values . 



11, A method according to Claim 9, wherein said 
10 evaluating step is to judge whether said result of the 
detection is correct or not, using a variance. 



12. An image processing method for judging whether 
an object area in an image includes an irradiation 
15 area, said method comprising: 

a secondary difference value acquisition step of 
acquiring secondary difference values from one- 
dimensional image data of said object area; 

an irradiation end extraction step of extracting a 
20 coordinate of an end point of said irradiation area 

from the secondary difference values acquired in said 
secondary difference value acquisition step; 

a comparison step of comparing the coordinate 
extracted in said irradiation end extraction step with 
25 a coordinate of an end point of the irradiation area 

included in said image, said coordinate being obtained 
preliminarily; and 
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a judgment step of Judging whether said object 
area includes the irradiation area, based on the result 
of the comparison in said comparison step. 



5 13. A method according to Claim 12, wherein said 

Judgment step comprises a step of Judging that said 
object area does not include the irradiation area, if 
the coordinates are close to each other, or otherwise 
judging that said object area includes the irradiation 



14. A method according to Claim 12, further 
comprising an accumulated image data production step of 
producing projection of said object area as said one- 

15 dimensional image data, 

wherein said irradiation end extraction step 
comprises a step of carrying out processing for the 
one-dimensional image data obtained in said accumulated 
image data production step. 

20 

15. An image processing method for judging whether 
an object area in an image includes an irradiation 
area, said image processing method comprising: 

a coordinate indication step of providing an 
25 indication of a plurality of rows for which one- 
dimensional image data is to be extracted from said 
object area; 
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a secondary difference value acquisition step of 
acquiring secondary difference values from the one- 
dimensional image data of said object area according to 
the indication in said coordinate indication step; 
5 an irradiation end extraction step of extracting 

coordinates of end points of said irradiation area from 
the secondary difference values acquired in said 
secondary difference value acquisition step; 

a storage step of successively storing the 
10 coordinates extracted in said irradiation end 
extraction step; 

an average acquisition step of acquiring an 
average of the plural coordinates stored in said 
storage step; 

15 a comparison step of comparing the average of the 

coordinates obtained in said average step with a 
coordinate of an end point of the irradiation area 
included in said image, said coordinate being obtained 
preliminarily; and 

20 a first judgment step of judging whether said 

object area includes the irradiation area, based on the 
result of the comparison in said comparison step. 

16. A method according to Claim 15, wherein said 
25 first judgment step comprises a step of judging that 

said object area does not include the irradiation area, 
if the coordinates are close to each other, or 
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otherwise judging that said object area includes the 
irradiation area. 

17. A method according to Claim 15, further 

5 comprising a second judgment step which is carried out 
based on the result of the judgment in said first 
j udgment step , 

wherein said second judgement step comprises a 
variance acquisition step of acquiring a variance of 

10 the coordinates stored in said storage step, a 

comparison step of comparing the variance obtained in 
the variance acquisition step with a predetermined 
value, and a judgment step of judging whether said 
object area includes the irradiation area, based on the 

15 result of the comparison in the comparison step. 

18. A method according to Claim 16, wherein said 
second judgment step comprises a step of carrying out 
each of the steps when said first judgment step results 

20 in judging that said object area does not include the 
irradiation area. 

19. A method according to Claim 15, wherein said 
irradiation end extraction step comprises a step of 

25 carrying out said extraction of coordinate, based on 
the sign, either positive or negative, of a primary 
difference value of said one-dimensional image data. 
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20. A method according to Claim 15, wherein said 
secondary difference value acquisition step comprises a 
step of carrying out said acquisition of the secondary 
difference values from said one-dimensional image data 

5 subjected to smoothing. 

21. An image processing method comprising a 
coordinate indication step of indicating a row for 
calculation of characteristic quantities of two- 

10 dimensional image data, a characteristic quantity 

calculation step of calculating said characteristic 
quantities from data of the row indicated in said 
coordinate indication step, an end point extraction 
step of extracting an end point of an object area from 

15 the characteristic quantities calculated in said 

characteristic quantity calculation step, an end point 
storage step of storing coordinates of end points 
extracted in said end point extraction step, a rotation 
angle indication step of indicating an angle of a 

20 rotation axis onto which the end points stored in said 
end point storage step are projected, an accumulated 
quantity calculation step of calculating projection of 
the end points stored in said end point storage step 
onto said rotation axis of the angle indicated in said 

25 rotation angle indication step, an accumulated quantity 
storage step of storing a projection quantity onto said 
rotation axis, calculated in said accumulated quantity 
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calculation step, and a rotation angle judgment step of 
judging a rotation angle of the object area from said 
projection quantity stored in said accumulated quantity 
storage step. 

22. A method according to Claim 21, wherein a 
start point of said rotation axis onto which the end 
points stored in said end point storage step are 
projected is placed at a barycenter of image data not 
less than a fixed density value. 

23. An image processing apparatus comprising: 
means for determining a plurality of areas 

arranged in a predetermined direction on an image and 
each having a predetermined shape; 

means for calculating a secondary difference value 
of density values representing the respective areas in 
said plurality of areas; and 

means for judging one end point of an irradiation 
area from said secondary difference values calculated 
by said calculating means. 

24. An image processing apparatus for extracting 
an irradiation area in an input image, said image 
processing apparatus comprising: 

means for detecting an irradiation end, based on a 
density distribution in each area, for a plurality of 
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areas in a desired direction in said image; and 

means for evaluating the result of said detection, 
based on the result of irradiation ends detected for 
each of said plurality of areas. 

5 

25. An image processing apparatus comprising 
coordinate indication means for indicating a row for 
calculation of characteristic quantities of two- 
dimensional image data, characteristic quantity 

10 calculation means for calculating said characteristic 
quantities from data of the row indicated by said 
coordinate indication means, end point extraction means 
for extracting an end point of an object area from the 
characteristic quantities calculated by said 

15 characteristic quantity calculation means, end point 
storage means for storing coordinates of end points 
extracted by said end point extraction means, rotation 
angle indication means for indicating an angle of a 
rotation axis onto which the end points stored in said 

20 end point storage means are projected, accumulated 

quantity calculation means for calculating projection 
of the end points stored in said end point storage 
means onto said rotation axis of the angle indicated by 
said rotation angle indication means, accumulated 

25 quantity storage means for storing a projection. 

quantity onto said rotation axis, calculated by said 
accumulated quantity calculation means, and rotation 
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angle judgment means for judging a rotation angle of 
the object area from said projection quantity stored in 
said accumulated quantity storage means. 

5 26. A storage medium storing a program for 

carrying out an image processing routine comprising: 

a step of determining a plurality of areas 
arranged in a predetermined direction on an image and 
each having a predetermined shape; 
10 a step of calculating a secondary difference value 

of density values representing the respective areas in 
said plurality of areas; and 

a step of judging one end point of an irradiation 
area from said secondary difference values calculated 
15 in said calculating step. 

27. A storage medium storing a program for 
carrying out an image processing routine for extracting 
an irradiation area in an input image, said image 
20 processing routine comprising: 

a step of detecting an irradiation end, based on a 
density distribution in each area, for a plurality of 
areas in a desired direction in said image; and 

a step of evaluating the result of said detection, 
25 based on the result of irradiation ends detected for 
each of said plurality of areas. 
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28 . A storage medium storing a program for 
carrying out an image processing routine comprising a 
coordinate indication step of indicating a row for 
calculation of characteristic quantities of,two- 
5 dimensional image data, a characteristic quantity 

calculation step of calculating said characteristic 
quantities from data of the row indicated in said 
coordinate indication step, an end point extraction 
step of extracting an end point of an object area from 

10 the characteristic quantity calculated in said 

characteristic quantities calculation step, an end 
point storage step of storing coordinates of end points 
extracted in said end point extraction step, a rotation 
angle indication step of indicating an angle of a 

15 rotation axis onto which the end points stored in said 
end point storage step are projected, an accumulated 
quantity calculation step of calculating projection of 
the end points stored in said end point storage step 
onto said rotation axis of the angle indicated in said 

20 rotation angle indication step, an accumulated quantity 
storage step of storing a projection quantity onto said 
rotation axis, calculated in said accumulated quantity 
calculation step, and a rotation angle judgment step of 
judging a rotation angle of the object area from said 

25 projection quantity stored in said accumulated quantity 
storage step. 
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ABSTRACT OF THE DISCLOSURE 

An image processing method is intended to extract 
an irradiation area in an image easily and accurately. 

A calculation unit 202 calculates secondary 
difference values of density values representing 
respective areas in three areas of a rectangular shape 
arranged in parallel on an image, determined by a 
calculation area input unit 200 and a calculation area 
determination unit 201, between adjacent rectangular 
areas and the secondary difference values thus 
calculated are stored in a memory unit 203 . Then a 
judgment unit 204 Judges one end point of an 
irradiation area from the secondary difference values 
thus stored and an irradiation area determination unit 
205 determines the irradiation area, based on a 
plurality of end points thus judged. 
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